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ABSTRACT 


A technique is presented for calculating the capacity of 
a lake for development based on quantifiable relationships between 
nutrient inputs and water quality parameters reflecting lake trophic 
status. From the land use and geological formations prevalent in a 
lake's drainage basin, the phosphorus exported to the lake in run- 
off water can be calculated, which, when combined with the input 
directly to the lake's surface in precipitation and dry fallout, 
gives a measure of the natural total phosphorus load. From the 
population around the lake, the maximum artificial phosphorus load 
to the lake can be calculated and can be modified according to 
sewage disposal facilities used, if such areeracian is warranted. 
The sum of the natural and artificial loads i.e. the total load to 
the lake can be combined with a measure of the lake's morphometry 
expressed as the mean depth, the lake's water budget expressed as 
the lake's flushing rate, and the phosphorus retention coefficient 
of the lake, a parameter dependent on both the lake's morphometry 
and water budget, to give a prediction of the springtime total 
phosphorus concentration in the lake. The latter prediction is based 
on a model reported by Vollenweider (1969). Long-term average runoff 
per unit of land area, precipitation and lake evaporation data for 
Ontario provide a means of calculating the necessary water budget 
parameters without expensive and time-consuming field measurements. 
The predicted spring total phosphorus concentration can be used to 
predict the average chlorophyll a concentration in the lake in summer, 
and this, in turn, can be used to estimate the Secchi disc transparency. 


Thus, the effects of an increase in development on a lake's water quality 
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can be predicted. Conversely, by setting limits for the "permissible" 
summer average chlorophyll a concentration or Secchi disc transparency, 
the "permissible" total phosphorus concentration at spring overturn 

can be calculated. This can be translated into "permissible" artificial 
Toad, which can then be expressed as total allowable development. This 
figure can be compared to the current quantity of development and 
recommendations made concerning the desirability of further development 


on the lake. 
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The relationship between Secchi disc and chlorophyll] a 
for thirty-five lakes in the Haliburton Highlands region 
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Mean annual precipitation a) Southern Ontario 

b) Northern Ontario. 


Mean annual lake evaporation (from Bruce and Weisman 1966). 
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GLOSSARY OF SYMBOLS 


Meaning 


lake area 

mean depth 

lake volume 

development of shoreline 

drainage basin area 

unit runoff 

precipitation 

lake evaporation 

total lake outflow volume 

areal water load = Q/A, 

flushing rate 

flushing time = 1/p 

phosphorus retention coefficient for 
lake 

phosphorus sedimentation rate 

loading 

"permissible" loading 

loading from virgin drainage area 

loading from precipitation 

natural loading = Le + Lop 

artificial loading 

total loading = Ly + Ly 

supply 

“permissible” supply 


supply from virgin drainage area 


Dimensions 
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Meaning 


supply from precipitation 
natural supply 

artificial supply 

total supply 

number of cottages 

number of permanent dwellings 
total capita-years per year 


additional "permissible" cottages 
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INTRODUCTION 


The demand for summer cottaging and other recreational 
facilities in Ontario is increasing rapidly, with about 10,000 new 
cottages under construction in each year (Dept. of Tourism and 
Information). Governmental personnel (planners, biologists, engineers ) 
having responsibilities of evaluating and subsequently approving or 
rejecting proposed developments are faced with the difficult task 
of effecting decisions without guidelines which quantitatively 
determine the effects or impacts of such development on the quality 
of the aquatic and terrestial environment. It must be stressed that 
the decisions made most assuredly have a direct effect on that quality. 
Suitable guidelines must be based on predictive ecological and social 
theories which, to date, have not been formulated in terms of practical 
management tools. It is the aim of this manual to provide a means of 
determining the capacity of a lake for cottage or permanent home 
development. Assuming that specific acceptable limits can be established 
for a number of parameters which reflect the quality of the aquatic 
environment, translation of these limits to permiss{ble cottage numbers 
or other development schemes can be achieved through application of 
techniques described in this manual. Conversely, if a decision is made 
to allow the development of a given number of cottages or large-scale 
condominiums, then the governmental agency responsible will be completely 
cognizant of the environmental ramifications of their decision. 

At this time, only the problem of lake trophic status is 
considered, with the equally important concerns of the effects of 
development on fisheries, wildlife and on human health not taken into 


account (there is, of course, considerable relationship between trophic 
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state and fisheries potential.) It should be stressed that the 
approach described herein is intended to serve as an interim 
solution until the results of the Lakeshore Capacity Study (which 
is currently being undertaken by personnel of Ministry of the 
Environment and Ministry of Natural Resources), have been obtained 
and translated into management principles. This latter study will 
take cognizance of health, fisheries and wildlife, as well as 
trophic status, in assessing the impact of development on watersheds. 

In understanding lake conditions, it is important to realize 
that the entire watershed and not just the lake or the lake and its 
shoreline, is the basic ecosystem unit. The terrestial and aquatic 
portions of any watershed are inherently linked with the gravitational 
movement of minerals in drainage waters flowing from the land to the 
water (ie. lakes, oceans), as the major terrestial-aquatic linkage 
(Likens and Bormann, 1974). It is important to note that the aquatic 
portion of a watershed is "downhill" from the terrestial portion; that 
is, dissolved and particulate material from the land are transported 
by geological processes to the water and eventually to the less accessible 
sites of the sediments in the aquatic portions of the ecosystem. There 
is, of course, some "uphill" movement of material by meteorological and 
biological transport, eg. release of No from a lake through denitrification 
and gas diffusion, followed by precipitation of the same to the land 
for the former process, and movement of material through the food chain 
for the latter (Figure 1). However, the net movement is always "downhill", 
ie. to the aquatic system, with the deficit made up by geological weathering 


of rocks. The implication of this in terms of a lake's capacity for 


development is that any alteration in a watershed ultimately, through 
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a_ complicated series of terrestial-aquatic linkages, effects 


the lake. Thus, any management approach based solely on the 


lake and its shoreline is both simplistic and inaccurate. 
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BACKGROUND FOR A TROPHIC STATUS INDEX 


Historical Development 


As stated in the Introduction, adequate watershed 
management can only be based on quantitative predictive 
ecological theories. Very few useful models or theories are 
available to the planner at the present time because ecology 
has been largely a descriptive science with too little attention 
given to predictability. Emphasis must now be directed to the 
development of models whose value will be judged on their predictive 
ability alone. Under these conditions, empirical models are as 
valid as realistic models provided they make testable predictions 
about the effects of various stresses on the environment. 

Early empirical approaches were directed almost exclusively 
to fisheries management in the sense of potential yields or production 
rather than to lake trophic state. Rawson (1952, 1955) related the 
25-year avérage commerical fish production for a number of large 
Canadian lakes to mean depth - an index of the lake's morphometry 
(Figure 2), a parameter recognized by Thienemann (1927) many years 
before as one of the most important factors controlling lake productivity. 
Moyle (1946) related fish production in Minnesota to total alkalinity 
and total phosphorus content; Rounsefell (1946) found an inverse 
relationship between lake surface area. Hayes and Anthony (1964) 
combined the above approaches, relating fish yield to lake area, mean 
depth and alkalinity. Other researchers, notably Larkin and Northcote 
(1958), Hrbagek (1969), and Reimers, Maciolek and Pister (1955), also 
investigated and reported similar relationships. 

In 1965, Ryder introduced what is now the most popular method 
of estimating potential fish production, the morphoedaphic index. This 


index, determined as the total dissolved solid concentration divided by 
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lake mean depth, was highly correlated to fish production for 23 
Canadian lakes. The morphoedaphic index has been further investigated 
by Jenkins (1968, 1970) and has been found to be suitable for a 

wider range of lakes than the original selection. The most 
satisfactory aspect of this index is that a prediction can be made 

for potential fish yield from only two simple parameters, the total 
dissolved solids content and mean depth of the lake. However, this 
approach is not suitable for determining how development will effect 

a potential yield. 

Although considerable attempts to manage fisheries by the 
empirical approach have been tried, very little has been accomplished 
in terms of managing development by a similar method. A notable 
exception was Seppanen's (1972) proposal to determine a lake's summer 


cottaging capacity. The author suggested that a suitable formula for 


the recommended number of shoreline cottages was wl where A is 

the lake area in hectares, Dj is the development of shoreline equal 

to the shorelength divided by the circumference of a circle of the 

same area, and 10 is a figure representing the minimum area of lake 
surface (in hectares) needed per cottage for a round lake with no 
islands. While providing a formula that takes into account both lake 
area and shoreline length, many other important factors are neglected 
making the technique quite limited for its intended purpose. Presently, 
in Ontario, governmental approval or rejection of a resort development 
is determined, in part, through application of individual components of 
the Lake Alert Study (Hough, Stansbury and Assoc. Ltd. 1972). The 
simplest and most frequently used procedure, the BOAT/LIMIT system, 

is based on a requirement of 3.3 acres of "usable" lake surface (usable 


in terms of power boating) per cottage. A much more sophisticated 
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system is also outlined, the SCORE/MAP system. Under this scheme 
constraints are mapped for the shoreline based on water quality, 

fish and wildlife, vegetation, soil, scenic elements, boating, 
access, private lands, and public outdoor recreation. However, 

the relative importance of each of these factors must be arbitrarily 
decided and a weighting factor employed. The constraints are summed 
and an arbitrary cutoff point is determined. Using a minimum lot 
size guideline, the acceptable shoreline is divided and the minimum 
of either this number or the boat limit result is termed the acceptable 
total development. Considering the water quality aspects of Lake 
Alert, six parameters (mean depth, oxygen distribution, chlorophyll a 
concentration, Secchi disc, iron:phosphorus ratio, and morphoedaphic 
index) were measured; the resultant values were scored between 0 and 
10. The sum (0-60) was then ranked along with other constraining 
parameters. The shortcoming of this approach is that no predictive 
relationship between development and potential water quality (or any 
other parameter) can be made, ie. a calculation that a given lake has 
a capacity for 200 additional cottages does not include any prediction 
of how those cottages will effect water quality or any of the other 
constraining parameters. The Lake Alert scheme is, in fact, based 

on qualitative rather than quantitative evaluation. 

An entirely different and far-reaching, although somewhat 
indirect, approach has recently been used in a lake management sense. 
In 1968, Vollenweider published in his classic text on eutrophication, 
a plot (Figure 3) of total annual phosphorus loading (the amount of 
material added per unit surface area of lake per annum) vs. mean 
depth (z). Vollenweider found that bands could arbitrarily be drawn 


separating the lakes into the three standard lake types in terms of 
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degree of eutrophy: oligotrophic, mesotrophic and eutrophic. 
The lower line separating oligotrophic from mesotrophic lakes 
was termed the "permissible loading" since it represented the 
upper loading level as a function of mean depth that could be 
permitted without having the lakes revert beyond the oligotrophic 
state, while the upper line, the "critical loading", represented | 
the level above which a lake could be characterized as eutrophic. 
A similar relationship was derived for nitrogen loading assuming 
a tolerable N/P ratio of 15:1 (by weight). 

Faith in this simple model stemmed from the fact that 
for lakes with phosphorus loading data available, the predicted 
status in most cases matched the observed status as exemplified 
by the standard (but arbitrary) criteria of lake condition: 
transparency (Secchi disc depth), chlorophyl1 a concentration, 
oxygen concentration in the hypolimnion, frequency of algal blooms, 
etc. This simple relationship has been widely accepted both as a 
guide to the degree of eutrophy of lakes and, more important, as 
a guide to the "permissible" and "dangerous" loading levels of 
phosphorus to lakes (Schindler and Nighswander 1970; Schindler 
et alae Gachter and Furrer 1972; Patalas 1972; Patalas and 
Salki 1973; Ahl 1972; Stumm and Stumm-Zollinger 1972; International 
Joint Commission 1969). Perhaps the most important use of this 
relationship was to provide a rationale for the removal of phosphorus 
from sewage treatment effluent in the Laurentian Great Lakes. For 
example, Lake Ontario, in the mesotrophic category according to its 
estimated total phosphorus loading for 1967, would, in 1986, fit into 
the oligotrophic range of the graph if 95% of the phosphorus was 
removed from all municipal and industrial wastes. However, without 
phosphorus control, Lake Ontario would, according to Figure 3 be well 


into the eutrophic range in 1986. By a similar argument, Lake Erie, 
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well in the eutrophic range in 1967, would be improved almost 
to the mesotrophic category in 1986 with phosphorus control or 
would be even more eutrophic without control. In much the same 
way, if a lake's phosphorus loading is measured or calculated, 
the additional load required to maintain a given trophic status 
level can be calculated and loosely interpreted in terms of 
capita-years; this calculation in turn may be translated into 
cottages or additional capacity for development. 

The first reported case of lakes which clearly did not 
fit Vollenweider's phosphorus loading - mean depth characterization 
was described by Dillon (1974a,b). In a study of the total phosphorus 
budget of 19 lakes in southern Ontario, it was found that a number 
of lakes had very high phosphorus loadings (5 lakes with L>1 g m-2 
yr7l), According to Vollenweider's relationship, many of these lakes 
should have been eutrophic, yet in all cases the chlorophyll a 
concentrations were low (summer averages between 0.8 and 2.7 mg m?), 
Secchi disc depths were high (summer averages, 3.5 - 8.7 m) and 
significant hypolimnetic oxygen deficits were not formed. The reason 
for this discrepancy was undoubtably that many of these lakes had 
very high flushing rates resulting from the large drainage area/lake 
area ratios. 

In 1973, Vollenweider modified his relationship to take 
such lakes into account. This was accomplished by incorporating the 
mean residence time of the water in the lake (t,). By plotting L vs. 
Zs (Figure 4), a more realistic representation was achieved. The 
parameter Z/ty is equivalent to the areal water loading (ie. the 
height of the water load over the lake's area that is supplied in one 
year). The L vs. ra model has been used for predicting changes in 


trophic status of several Muskoka Lakes (Michalski, Johnson and Veal 
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1973) that would result from changes brought about by the 
reduction of phosphorus in sewage effluent (Figure 5). As 
effected by Michalski et al (1973), the model can be employed 
by environmental managers in fie Same manner as the L vs. z 
plot: the maximum additional loading to a lake that will not 
surpass the "permissible" loading can be calculated and inter- 
preted in terms of cottage development. Although a significant 
improvement over the L vs. z plot, this method still lacks 


quantitative predictability in terms of water quality parameters. 
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Basis for a New Approach 


Any approach to predict water quality from a trophic 
status point of view must take into account the importance of 
the element phosphorus. That phosphorus is the nutrient most 
frequently controlling production and therefore trophic status 
in north temperate lakes has been demonstrated on numerous 
occasions (Schindler et al. 1971, Schindler 1974; Fuhs et al. 
1972). As early as 1947, Sawyer (1947) recognized that phosphorus 
concentration in the lake water was the factor controlling 
eutrophication. Although Vollenweider's L vs. Z plot relates 
loading rather than concentration to trophic state, his improved 
second relationship, L vs. Z/ es relates phosphorus concentration, 
and not loading, to trophic state. This may seem at first to be 
a contradiction but the rationale is as follows: the lines separating 
lakes into distinct trophic types have the dimensions of U/ Za or 


g m-2 yr! /m yr7l, jie. g m3, 


Thus, the lines are independent of 
time and have units of concentration (the lines of the L vs. z plot 
have units of g m73 yr, ie. of volumetric rather than areal loading). 
Since the parameters with which one subjectively evaluates a lake 
(chlorophyll a, oxygen in the hypolimnion) are also expressed in 
terms of concentration, the above interpretation is most reasonable. 
Consequently, predictions relating the impact of development on the 
phosphorus concentration of a lake and subsequently on parameters 
describing the trophic state are central to a predictive management 
scheme. | 

The overall approach employed in this manual is shown in 


Figure 6. From consideration of the geology and land use of a lake's 


drainage basin, it is possible to estimate the total phosphorus 
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exported or washed out per unit area of watershed, which, combined 

with the drainage area provides an estimate of the total phosphorus 
Supplied to the lake from the land. Addition of the input of 

phosphorus in precipitation directly falling on the lake allows 
calculation of the natural phosphorus load to the lake. Development 
existing on the lake is then measured (aerial survey or field counts) 
and the phosphorus loading from artificial sources calculated. The 
total loading, natural plus artificial, may then be combined with the 
lake's morphometry and water budget to predict a phosphorus concentration 
which is subsequently related to the average summer chlorophyll a 
concentration. From this latter calculation, one may determine the 

mean ice-free Secchi disc visibility. This manual will most often be 
used in the opposite sense. For example, maximum acceptable average 
summer chlorophyll] a concentration (or minimum Secchi disc reading) 

will be established by the appropriate governmental agency. From these 
limits, decision-making personnel will be able to calculate a maximum 
permissible phosphorus concentration, which can be interpreted in terms 
of a maximum "permissible" total phosphorus load, because the lake 
morphometry and water budget are essentially fixed. That is, although 
there are year-to-year variations in the lake's water budget, the long- 
term average budget is acceptable. Finally, the maximum permissible 
artificial phosphorus loading can be estimated and expressed as the 
maximum allowable development (ie. numbers of cottages, etc.). Although 
hypolimnetic oxygen deficits have been related to phosphorus loading for 
the Great Lakes (Gilbertson, Dobson and Lee 1972), a direct relationship 
applicable to a wide variety of lakes and especially to lakes of the 
Precambrian Shield is not yet available. Mortimer (1941, 1942) suggested 
limits of 0.025 mg cm™¢ day”! for eutrophic lakes, but the quantitative 
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link to nutrients is, as yet, unformulated. 

In addition, the potential fish yield of the lake to be 
considered for development can be approximated by Ryder's morphoe- 
daphic index (TDS + Zz), but the effect of the development on future 
fish yields cannot. The change in TDS resulting from a change in 
the nutrient budget brought about by development cannot be predicted. 

In the following section the theory and method for each 
step of the scheme are described in detail, and subsequently a 
stepwise procedure for calculating the lake's Capacity is outlined 


for use by planners, managers, etc. 
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THEORY FOR EACH STEP OF THE MANAGEMENT SCHEME 


Calculation of the Phosphorus Loading to the Lake 


1. Natural Phosphorus Load (Ly) 


(a) From Land (Le) 


The geological formations of southern Ontario can, as a 
first approximation, be classified as Precambrian igneous rock of 
plutonic origin (Canadian Shield) or sedimentary rock. The former 
is typically composed of granites, gneisses, pegmatites, syemite, 
migmatites, diorite, gabbro, hornblendite, amphibolite and pyroxenites, 
the latter of limestone, dolomite, shale and basal clastics. A very 
large proportion of Ontario's recreational lakes are situated on the 
Shield, with only the southern portion of the Kawartha - Trent 
system, Lake Simcoe and a few smaller lakes on sedimentary material. 
In addition, the watersheds of most of the Shield lakes are entirely 
or almost entirely forested, with the remainder being either marsh- 
land or pastureland which, although used for agricultural purposes, 
is not chemically fertilized. Even south of the Shield on the 
sedimentary bedrock, there are few areas where intensive agriculture 
(ie. chemical fertilization) is practised in prime recreational land; 
Some areas around Lake Simcoe are a notable exception. Dillon and 
Kirchner (1974) have developed a phosphorus export scheme (Table 1) 
which is based on classification of geology as either "igneous" or 
"sedimentary" and land use as "forest" or "forest plus pasture"; 
the latter category implies that 15% or more of the watershed is 
cleared but unfertilized land. This export scheme is applicable 
to almost the entire recreational lake area of Ontario. The results 
of Dillon and Kirchner are based on a combination of a Study of 
the phosphorus export of 34 watersheds in southern Ontario and 


all additional phosphorus export studies reported in the literature 
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Table 1: Ranges and mean values for export of total phosphorus 
from 43 watersheds. Results in mgm = yr 


Geological Classification 


Igneous Sedimentary 
Land Use 
Forest 
Range 0.7 - 8.8 6.7 - 18.3 
Mean 4.7 NO ea? 


Forest + Pasture 
Range 5.9 - 16.0 Hitele 37 40 
Mean 10.2 yA TE: 
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where watersheds fall into the above mentioned categories. It is 
important to note that the natural export from sedimentary materials 
(off the Shield) is almost exactly double that from igneous bedrock 
(on the Shield); therefore, lakes such as the Kawarthas have a 

higher natural loading than Precambrian Jakes. In addition, a change 
from a land use of forest to forest plus pasture doubles the export 
within a geological classification; for example, from 4.7 to 10.2 


mg m2 yr! for igneous watersheds and from 11.7 and 23.3 mg m™2 i) le 


in sedimentary watersheds. If a development significantly alters the 
amount of cleared land in a watershed, then an appropriate change in 
the export value used must be made. 


To calculate the natural phosphorus loading to a lake from 


] 


its drainage area, one must know the area of the watershed of each 
tributary to the lake, and be able to classify each as to geology 

and land use so that the export coefficient can be determined. Should 
the location be one of the few where intensive farming is undertaken, 
the reader should consult Vollenweider (1968), Dillon and Kirchner 
(1974) and Loehr (1974). The total amount of phosphorus supplied to 
the lake from the land is therefore calculated as the sum of the area 


of each drainage basin times its phosphorus export coefficient: 


Je = 2 (Ay° Es) 


The areal loading of phosphorus supplied to the lake from the land 
Les is equal to (Z Ag Es )/A, and is equivalent to the supply divided 
by the lake area. 

A complicating factor arises if any tributary or watershed 
of the lake in question has an additional lake or lakes in its course. 


These lakes undoubtably act as traps for phosphorus and other nutrients, 


decreasing the actual amount of material transported from the drainage 
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area to the lake in question. A means of accounting for this is 


described in a later section. 
(b) From Precipitation (Lop) 


Phosphorus input via precipitation (wet and dry fallout) 
has been virtually ignored until recent years. Studies by Schindler 
and Nighswander (1970), Armstrong and Schindler (1971), Barica and 
Armstrong (1971), Dillon (1974a), Dillon and Rigler (1974a), Kluesener 
(1972) and Shannon and Brezonik (1972) have demonstrated that, for 
many lakes, precipitation can be a major source of nutrients. From 


consideration of the above-mentioned studies, which found loadings 


ranging from 27 - 102 mg m2 of lake surface yr’! a value of 75 mg m7? 


yr! is recommended as applicable to southern Ontario lakes. Most of 
the lower values did not include a measure of phosphorus in dry fallout 
and the highest value, 102 mg m7? ae was measured in an area where 
much of the land was used for agriculture. Furthermore, the figure of 
77 mg m-e yr! found by Dillon (1974a) was measured in the Haliburton 
Highlands area of the province, an area of prime concern for recreational 
development. Thus, a figure of 75 mg m7? yr”! is reasonable, although 
additional research is obviously needed. This topic is reviewed by 


Chapin and Uttormark (1973). 


2 


With Lop, = 75 mg m yo the natural phosphorus supply to 


a lake (Jy) is given as: 


UN i Je ta) 


N PR 


= (x(Ay*E,) + 75-A,)/10° (kg yr!) (2) 


or the loading (Ly) 
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patente (mg m- yr!) 


Late + 75 
N 
Ao 


with A, and Aq in ne. 
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2. Artificial Phosphorus Load (Ly) 


The calculation of the phosphorus supplied to a lake by 

the population in its drainage basin is a difficult task and must 

necessarily be based on a supply per capita-year figure with 

several assumptions inherent in this method. Although numerous 

figures are quoted in the literature (eg. see Vollenweider 1968 or 

Dillon 1974c for a summary) great care must be taken in selecting 

the appropriate value. The following points must be remembered: 

a) Values measured in studies such as that of Johnson and Owen in 
1971-(185eke capita”! yy and used by Michalski, Johnson and 
Veal (1973) are no longer applicable because of legislation to 
reduce the phosphorus content of laundry detergents. High 
phosphorus content dishwashing detergents, however, remain 
legal and a survey of the Muskoka lakes found that 30% of the 
cottages in the area employed automatic dishwashers. 

b) The waste disposal technique for most developments that are 
already established or will be established in the near future 
is the conventional septic tank - tile field system. The 
efficiency of this treatment as far as phosphorus removal is 
concerned is dependent on the type and depth of soil surrounding 
the tile bed and between the tile bed and the lake. In Precambrian 


areas, typically having very shallow, coarse-textured sandy or 


muck soils there is no satisfactory evidence which indicates that 
phosphorus is retained in the soils. Therefore, it must be assumed 


that all phosphorus discharged to soils of a tile bed area eventually 
gains access to the lake. In sedimentary areas, septic tank - tile 
field systems located in sand, gravel or muck areas are likely to 


be as ineffective as far as phosphorus retention is concerned as 
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those on systems located on the Shield. Lakes surrounded by 
clay or clay-loam soil, however, will be provided with some 
measure of protection. 

c) Additional sources of phosphorus, (eg. fertilizer applied to 
lawns in cottage areas), are impossible to evaluate quantitatively 
without a lake by lake survey. 

With the above factors in mind, a reasonable phosphorus 

Supply per capita-year can be calculated. According to the data of 

Bucksteeg (described in Vollenweider 1968), the annual per capita 

amount of phosphorus uptake in food in Germany in 1960 was 0.55 kg.; a 

Similar amount is assumed to be excreted. In addition, the yearly per 

capita phosphorus supplied as domestic sewage (excrement + household 

wastes) in 13 studies in North America and Europe averaged 0.80 kg. 


be the studies having been carried out before use of high phosphate 


year. 
detergents was common. Thus, the two results are in basic agreement, 
and with food consumption in North America being greater than that in 
Europe and there being continued usage of high phosphate dishwashing 
detergents in North America, the higher figure of 0.80 kg phosphorus 
capita” | year”! is a good estimate of the phosphorus supplied from 
domestic sources. 


For all lakes situated on the Precambrian Shield the artificial 


supply of phosphorus can therefore be calculated as: 


= 70,6 .kg T capita-year 
Ja (capita-year Amit, Collages: x cottage-year ) (3) 


where N = the number of cottages or dwelling units 
T = the average number of capita years spent at each cottage in 
each year. 


The same formula can be used for areas south of the Shield that consist 
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of sand, gravel or muck soils. For clay or clay-loam soils in the 
sedimentary area, a supply of 0.4 kg capita”! year”! which takes 
into account the retention capacity of the soil can be used. 

The number of cottages or permanent dwellings (N) situated 
on a lake can be obtained by aerial photography (a good method is 
outlined in the Lake Alert Phase 2 report) or by field count. The 
time (in capita-years) spent at each dwelling per year (T) can be 
calculated using data gathered by the Department of Tourism and 
Information (1971). For cottages in the Haliburton region, 230 
capita days per year, (or 0.63 capita-years per year), and in the 
Kawarthas, 253 capita days per year (or 0.69 capita-years per year) 
are spent at each cottage (based on a survey figure of 4.3 people 
per cottage). For permanent dwellings, the appropriate value for 
T is.1,570 capita-days per year (or 4.3 capita-years per year.) An 
actual field survey for a particular lake would undoubtably improve 
these figures; nevertheless, they will serve as adequate, generally 
applicable estimates. 

In cases where holding tanks are employed and all wastes 
are removed to a treatment plant outside the lake's watershed, the 
supply of phosphorus from such dwellings should not be included in 
the calculations. 

Recent investigations involving the use of an iron-clay fill 
around the tile bed have demonstrated a potential phosphorus removal 
mechanism. In cases where such modifications to conventional tile bed 
systems are used and known to be in working order, a per capita 
phosphorus supply figure of 0.16 kg is suggested (80% removal). A 
factor, (1-R.), can be included in the model to allow for alterations 


in the phosphorus removal capacity of the system. 
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The total supply to the lake can be calculated as the 
sum of the natural and artificial supplies; 
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where R. is the retention capacity of the disposal system and will 
be equal to 0 for most Shield areas. The loading is, of course, 


equal to the supply per unit of lake area, 
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Prediction of the Phosphorus Concentration in a Lake 
1. Background 


The description of the theory and derivation of nutrient 
budget models would be both voluminous and mathematically difficult 
for most readers of this manual; for those interested, a review of 
these models published to the end of 1972 has been prepared by 
Dillon (1974d). Additional work by Vollenweider (1973) and Imboden 
(1973, 1974) is also highly pertinent to the topic. The model used 
in this manual derived from Vollenweider (1969) is employed because 
a) it is simply derived and can be used with a minimum of field 
measurements and b) it alone has been tested for lakes in southern 
Ontario and has been found to have good predictive capabilities. 

It can be assumed that the change in the concentration of 
phosphorus in a lake with time is equal to the supply added per unit 


volume minus the loss through sedimentation and the loss by outflow: 


CELE eee regs he opr te (5) 
dt V 


where [P] represents the phosphorus concentration (mg m?), J is the 
amount of phosphorus supplied per annum (tg yr7ly, V is the lake volume 
(m°), o is the sedimentation rate (yr!) and o is the flushing rate 
(yr) which is equal to the total volume of water outflowing per 

year (Q) divided by the lake volume (V). Since the volume of the lake 
is equal to the surface area (A) times the mean depth (z), equation (5) 


can be rewritten: 


aufhie 


dt 


eo tee LT (6) 


Net [rt 


where L = J/A, 
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The solution to this differential equation describing the concentration 


at a time t is: 


eee 1 = (1 - 2lote) EPI, ) g-(oto)t (7) 
z(o+p) L 


where Lae is the initial concentration. The Steady-state solution 


(t + ~) that we are interested in is simply: 


[p] = t 
Z(otp) (8) 


2. Use of the Equation 


Equation (8) can be used to predict the phosphorus concentration 
in a lake. The parameters required for this are the loading (L), mean 


depth (z), sedimentation rate (c) and flushing rate (p). 
(a) Loading 
Loading is calculated as described in the preceeding section. 
(b) Mean Depth 


Mean depth (as well as area and volume) for many lakes in 
Ontario can be obtained from the Ministry of Natural Resources Lake 
Inventory. If unavailable from this source, the lake will probably have 
to be echo-sounded, a contour map drawn, planimetered and the necessary 
morphometric parameters calculated. This is straightforward and is 


described in any standard limnological text (eg. Hutchinson 1957). 
(c) Flushing Rate 


The flushing rate (p) is calculated from the total outflow 


volume per year (Q) and the lake volume (V). Since this manual is 
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intended for use without fieldwork if possible, a method of calculation 

of Q is required. The long-term average annual areal runoff (r, in m yr!) 
has been mapped for Ontario south of the Shield by Coulson (1967) and 

for the entire Great Lakes Basin by Pentland (1968). Thus, virtually 

the entire area of recreational lakes in Ontario is mapped in terms of 

r, where r represents the differences between precipitation and 
evapotranspiration. The total long-term average inflow via surface 

runoff to a lake can therefore be calculated as Ayr. The total water 


balance can now be described by the equation: 


Q = Ayr baie (Pr-Ev) (9) 


where aval 8 represents the direct input of water to the lake by 
precipitation and Font Ey represents loss from the lakes by evaporation. 
The average value for Pr for a given area can be obtained from the 
Monthly Record or from isohyet maps (eg. Canada Land Inventory 1966) 

of Ontario. Maps of Ev (evaporation from lake surface) are given for 
Canada by Bruce and Weisman (1966). If Aa is large compared to Ay? 

then Q can be simply approximated as Bane: In the few studies undertaken 
in Ontario (eg. Schultz 1950), groundwater has made a negligible 
eonEribettan to the overal budget. The flushing rate, p, is then 


calculated as Q(m> yr” !y/v(m). 
(d) Sedimentation Rate 


Measurement of o, the sedimentation rate coefficient, for 
phosphorus in a lake is, at best, extremely difficult. Fortunately, an 
alternate parameter, the retention coefficient, R, can be shown to have 
a functional relationship to o. R is much more easily measured than o; 
(R is equal to the fraction of the loading that is not lost via the 


outflow) as well, the retention coefficient has been shown to be predictable 
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(Kirchner and Dillon, unpublished studies). 
The relationship derived between R and o by Dillon and 
Rigler (1974a) under the same assumptions as employed in deriving 


the original model is: 


Therefore, the equation (8) used to predict phosphorus concentration 


can be rewritten: 


[py] = L-U-R) 


Zp 
Kirchner and Dillon (1974) determined by multiple regression analysis 
that R was highly correlated with Q/A,, the areal water loading, usually 
written as qg. Their prediction of R is: 
R = 0.426 exp (-0.271 q,) + 0.574 exp (-0.00949 q,) (12) 

In some cases, the lake in question will have one or more lakes 
upstream which are sufficiently large to retain a significant amount of 
the total phosphorus exported from their respective portion of the 
watershed (see Figure 7). In the diagram the supply of phosphorus to 
Lake B is influenced by Lake A. This is taken into account by calculating 
the supply to Lake A, its retention coefficient (Ra), and multiplying Ja 


by (1-Ra) to give the fraction transferred to Lake B. 


3. The Concept of Response Time of a Lake 


As pointed out earlier change in the phosphorus loading to a 
Jake results in a change in the phosphorus concentration and a subsequent 


change in water quality. The change is, of course, not instantaneous but 
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Figure 7: 


Calculation of the natural supply of phosphorus to Lake B from 
its drainage area. Lake A acts as a trap for nutrient exported 
from area 3. A. does not include the area of Lake A. Additional 
input to Lake B’from precipitation and from artificial sources 
will be included in the final calculation of the supply. 
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rather is described in equation (7) and shown in Figure 8. 

The time required for a lake having an initial loading 
of Ly to respond to a change in loading to a new level Los (Lo can 
be greater or less than L,) can be described by the half-life of the 
change in concentration, ty (ie. by the time required for the lake's 
phosphorus concentration to move half-way from the original steady- 
State concentration to the final Steady-state concentration. Because 
equation (7) is exponential, two half-lives (2t, ) are required for the 
lake to reach 75% of its final concentration, 3° ty for 874%, etc. The 
half-life, ty. depends only on the rate coefficients o and OQ, representing 
the losses of phosphorus by sedimentation and outflow; it is independent 
of the loading level or the initial phosphorus concentration in the lake. 


It can be easily shown that 


t, = In 2/(otp) 


1 
2 


or for our use: 


» = {1-R) In 2 _ 0.69 (1-R) 
3 


B A (13) 


Therefore, one must consider the response time of a lake 
when predicting the effects of an increased loading (eg. resulting from 
development) or a decreased loading (eg. from improved sewage treatment 
facilities). It is suggested that 3 to 5 times the ty (ie. 87.5 - 96.9% 
of the way to the final steady-state concentration) be used as an 
indicator of the lake's response time. Characteristically, lakes with 
a rapid flushing rate will have short half-lives and therefore response 
times, while lakes that are very slowly flushed may take a long time to 
respond to a change in loading. For example, a lake with flushing rate 
of p = 6 yr"! (ie. the lake's volume is replaced 6 times per year by 
flushing) and retention rate of R = 0.3 has a ty of 0.69 (.7)/6 = 0.08 yr. 
Therefore, between 0.24 and 0.40 years are required for the lake to approach 


a new steady-state following a change in loading. 
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Relationship of the spring Phosphorus Concentration to the 


Summer Chlorophyll a Concentration and Water Transparency 
(Secchi Disc) 


Based on the work of Sakamoto (1966), Dillon and Rigler 
(1974b) developed a predictive relationship suitable for estimating 
the average summer chlorophy11 a concentration in lakes with spring 
N:P ratios > 12 (Figure 9). The data comprising this relationship 
were based on a combination of information presented by Sakamoto 
for 30 Japanese lakes, results for a wide variety of North American 
and European lakes reported by various authors, and information 
for 19 lakes in southern Ontario Studied by Dillon and Rigler. The 
equation for the prediction of the summer average chiorophyll a 


concentration from the spring phosphorus concentration is: 


10949 [ch] a] = 1.45 log, 9 [P] - 1.14 (14) 
with [ch] a] and [P] in mg m3 
North temperate lakes typically have NsP ratios of 15 - 40. 

The value of this relationship or "model" has been 
demonstrated in several cases. To date, the most comprehensive 
study of the effects of reduced nutrient load on a lake is that of 
Lake Washington (Edmondson 1969a). In 1952, ten sewage treatment 
plants discharged directly into the lake. Conditions in the lake 
(chlorophyl1] a concentration, Secchi disc depth, hypolimnetic oxygen 
deficit, phosphorus concentration, etc. ) gradually deteriorated 
(detailed by Edmondson 1961, 1966, 1968, 1969b, 1972; Anderson 1960) 
until 1963, the year following initial steps to divert sewage. 
Diversion was completed by mid-1968, thereby reducing the phosphorus 
loading by 55%, at which point improvements were already noticed. 


Mean chlorophyll] a concentration (July and August) which had been 


, ’ ; 
| : uJ @ i. i 


‘ _ 
a a 


| 


5 
gar gaitrg? Pe ih ‘afilen pea’ 


meee = 


? ye 
‘ey a 
20 oa wi ibs itiaes 5 1 Ug | eg es a 
—— - = _ 
€ 
= * : _ ' in 
shOGh UN | oe) Loetgor tana” 
; : ~ : ‘ BS Pee ol Ae Fe ‘ 
ae : : 
er | 
; : i 
"aloth Goeyieshet, (O51): atupsise tos > NO bSeBR, 
_ Ty 
4 tie » = iy | Woe t re "\ ¢ Se i ae } - Te Wh (Aaye 
4 : = : \= is _ vir 
ie 1S APE 25ABE OT MOTTSUIMSTAOD So PYNQSTOTET: aud SPs 739Vh riba 7 He 
! > a! : 
bids PAS. Dt? nd Oo ayunl 
‘ { 7 5 & 2% é \ q ‘i 4 
3 L ' ts t “ i { 
7 
G — 
; ” 1a a i oY 
| 
1" i y . { i ar »7 
; i y i 
= > t | + : 
! 4 é 
' - 
3 . : 
‘ ts ; '} 
rd aus oad | 
pd : - 
‘4d , : {) ‘ 
i ! t - , ‘ } i t Ne 
‘ 7h ri Py - } | ( \ ‘ I, * fio ( >| ‘ i - 
f | 7 : 
; } md 73 Ty oh » Bl yy ] oO) , a) : 
one Siti : 2th i “ans Gon ; ) 
A: (XO aT3 barG 1 fia orth rejo8 t em T 
PF ia'y ; P 
boo (BTteS (Tea } Loan sOT JENS ag 
' . 
a | F 
reer ¢ ae)! See ) 4 34 ~¥ ye a 
(OGG! aovraitA , 5 9eGl et Ri GET dae) det. soebianbd vcd 


- Sptiwne raw we 


Ps 


re ooh ad OMiwal fe) in 


35. 


100 


sar 
© 


JAPANESE LAKES 

a OTHER LAKES DESCRIBED 
IN THE LITERATURE 
HALIBURTON 


SUMMER CHLOROPHYLL a mg m7 


oa 
oO 


(SAKAMOTO, 1966) 


0-1 | | 
10 100 1000 


TOTAL PHOSPHORUS mgm? 


Figure 9: Relationship between average summer chlorophy1] a concentration and 
spring total phosphorus concentration (from Dillon and Rigler 1974b). 


MOTAUBIIAN -© 


- 
“Vy 
or 


fre 


BB: ‘norte 
ead 44 


cenwe 


about 1.5 mg m7? 


$ 


in 1950 and had increased steadily to about 
24 mg m° in 1964, was reduced to 3.5 mg m™> in 1969. The 
important point is that the changes in phosphorus and chlorophyll a 
very closely parallelled those predicted by the Sakamoto plot 
(Figure 10). 

In 1971, Little Otter Lake, a small Precambrian lake in 
Ontario, received a relatively large amount of a polyphosphate 
de-scaling agent (Michalski and Conroy 1973). The lake quickly 
developed classical symptons of eutrophy including a dense bloom 
of Anabaena limnetica G.M. Smith, low Secchi disc depths (0.1 - 
2.3 m), and high chlorophyll] a (12 - 46 mg m7?) and total phosphorus 
(12 - 110 mg m->) concentrations. Polyphosphate discharges were 
discontinued in the fall of 1971; phosphorus concentration in the 
following spring dropped to 11.5 mg m7? which, according to equation 
(14) should have resulted in an average chlorophy11 a concentration 


during the ice-free period of 1972 of 2.53 mg m->, 


The measured 
average concentration was 2.52 mg nm: an almost perfect fit (see 
Figure 10). 
| Additional phosphorus - chlorophyl]l a data have been supplied 
by Scheider and Rigler (unpublished studies) for 6 lakes in Algonquin 
Park, by Oglesby and Shafner for some of the Finger Lakes in New York 
State, and by Bachman and Jones for some Iowa lakes. These data are 
also shown in Figure 10 and it is apparent that in all cases the fit 
to the model is very good. 

It is expected that chlorophyl1 a concentration will be 
inversely related to the Secchi disc reading, (a measure of the water's 
transparency). Collection of data for a large number of lakes is shown 


in Figure 11. From the predicted chlorophyll a concentration, a 


prediction of the likely average Secchi disc reading is possible. Care 
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Figure 10: Additional spring phosphorus - summer chlorophyll data with line as 
derived from data shown in Figure 9. 
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SECCHI DISC (meters) 


O 5 10 15 20 25 
CHLOROPHYLL a (yg/!) 


The relationship between Secchi disc and chlorophyl] a for thirty- 
five lakes in the Haliburton Highlands region of Ontario. Values 
for each lake are based on means of values collected during the 
summer of 1973. Also, information from a number of other lakes is 


included as an indication of the relative status of lakes in the 
Study area. 
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must be taken in interpreting this relationship in the case of 
dystrophic (brown water) lakes, whose Secchi disc readings are 
lower than would be expected on the basis of chlorophyll a 


concentration alone. 
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Potential Fish Yield 


Ryder's (1965) morphoedaphic index (MEI), TDS/Z (total 
dissolved solids divided by mean depth) provides the best method 
of estimating the potential fish yield of the lake in question. 
The concentration of total dissolved solids in a lake fluctuates 
little over the year, so the time of measurement is not critical. 
The equation used to predict yield is: 


0.446 


F = 1.10 (TDS/Z) (15) 


where F is the fish yield in kg/ha, TDS is measured in ppm and z 
in meters. 

The limitations of this model should be stressed. Although 
an estimate of the potential sustainable fish yield is given, the 
effects of development on the fishery is unspecified. Over-develop- 
ment leading to excess nutrient input followed by a winter "fish 
kill" due to oxygen depletion under the ice or followed by greatly 
increased fish growths rates obviously results in a change in the 
fishery, but there is not necessarily a change in the MEI that 


would lead to prediction of either case. 
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STEPWISE PROCEDURE FOR CALCULATING THE DEVELOPMENT CAPACITY OF A LAKE 


Src it. 


Based on long-range plans for the lake, decide what the 


will be: 


Level 1: 


Level 2: 


Level 3: 


Level 4: 


maximum permissible average summer chlorophyll a concentration 


2 mg ms for lakes to be used primarily for body 
contact water recreation, and where it is desirable 
to maintain hypolimnetic concentrations of oxygen 


in excess of 5 mg ee 


The lake will be extremely 
clear with a mean Secchi disc visibility of >5 m 

and will be very unproductive. (Note - The Secchi 

disc visibility may be lower in brown water (dystrophic) 
lakes). 

5 mg m™>; for lakes to be used for water recreation 

but where the preservation of cold water fisheries 

is not imperative. The lake will be moderately 
productive and correspondingly less clear, with a 

mean Secchi disc visibility of 2 - 5m. 

10 mg m7: for lakes where body-contact recreation 

is of little importance, but emphasis is placed on 
fisheries (bass, walleye, pickerel, pike, maskinonge, 
bluegill, yellow perch). Hypolimnetic oxygen depletion 
will be a common occurrence. Secchi disc depths will 
be low (1-2 m), and there is a danger of winterkill 

of fish in shallow lakes. 


3, suitable only for warm water fisheries. 


25 mg m 
Secchi disc depth <1.5 m, hypolimnetic oxygen depletion 
beginning early in summer, considerable danger of 


winterkill of fish except in deep lakes. 
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The planning agency may pick any intermediate level should 
it so desire. 

STEP 2. From the chosen summer average chlorophy11 a concentration, 
calculate the permissible spring phosphorus concentration, 
[P] from: 
log fchiiaj ayl.45 logy, [P] - 1.14 


109, Lohivaereal . 14 
je. 10919 [P] = ——_— 


1.45 
eg. [chi a] [P] 
2 mg fies 9.9 mg m > 
5 18.5 
10 29.9 
Bs) 507.3 


STEP 3. Determine the lake surface area (A, in m), mean depth (z in m) 
and volume (V in m°) from the Ministry of Natural Resources 


be 


Lake Inventory (multiply acre-feet by 1233.6 to get m 
such data are not available, the lake must be properly sounded 
(preferably with an echo-sounder) and a contour map drawn. The 
area (A) is obtained by planimetry from an aerial photo of 

Known scale. 

STEP 4. Outline the lake's drainage area on a 1:50,000 scale topographic 
map and calculate the area (A, in m°) by planimetry (example 
shown in Figure 12). The following guidelines may help delineate 
the drainage basin boundary: 

1) boundaries are generally drawn by following the high points 
of land 
2) assume that water flows downhill to the nearest contour line 


3) assume that groundwater follows the drainage pattern of surface 


water. 
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Example of a lake's drainage area outlined from a 


topographic map. 
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STEP 5. From Figure 13 taken from Pentland (1968) determine the 


total amount unit runoff (r) in cfs min? and convert to 
m° yr”! m* or m yr’! by multiplying by 0.345. 


STEP 6. Calculate Q, the total outflow volume as: 


3 


A ere in yr!) 


d 
and calculate the flushing rate (po) as Q/V or 


CAA oir) /¥ (yr!) 
Gy Ay isy< 10 Ay? determine the mean annual precipitation 
(Pr) from Figure 14 and the mean annual lake evaporation 
(Ev) from Figure 15, convert to m yr”! by multiplying by 
0.0254 and calculate Q: 


Q=A,: a sa (Pr-Ev) 


d 
AL, on feet Py By) 
“ 
Ae eg ae (yr ) 
V 
STEP 7. Calculate the areal water load (q.) as 
V/A, (m yr!) 


STEP 8. Calculate the retention coefficient (R) as 
R = 0.426 exp (-0.271 q,) + 0.574 exp (-0.00949 q,) 
STEP 9. Calculate the response time of the lake to a change in 
phosphorus loading: 


Response time = (3 > 5) t 


nw 


= (3 + 5) (1-R) In 2 
p 


(yr) 


This will provide an indication of the time required for a 
lake to "respond" to development and will give an idea of 
when follow-up studies (if any) should be carried out. 
Conversely, the response time can assist in the interpretation 
of present water quality. For example, a lake with 250 new 
cottages may appear to be in good condition, but, if one 


can calculate that it has a response time of 6 - 10 years, 
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Figure 14b: Mean annual precipitation - Northern Ontario. 
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SLEEP “1:03 


SEP ds 


SrEP old, 


eae 


then caution is necessary before additional development 


is allowed. 
Calculate the permissible phosphorus load oor to the lake: 
ee Py ah ez. oe -2. =] 
Laie oe Pa ak (mg m~ yr “) 
det iSsi 
and the permissible supply oer 
J a ee) 2 
perm BREVIS 0 (kg yr -) 


10® 
Divide the watershed of the lake (on the topographic map) 
into subunits for all inflows and determine the area for 
each one (Anaad From the most recent aerial photos available 
determine if > 15% of the area of each watershed is either 
cleared land or marsh. Determine whether the watershed 
(not just the lake) is on Precambrian igneous rock or on 


aD 


sedimentary rock, and estimate the export (E mg m 
for each subwatershed from Table 1. Calculate the total 
supply of phosphorus from the land to the lake: 


2 : 6 -| 


and the load 


an 5 ay Nv at 
Lea - Aa E)/A, (mgm yr) 


Take consideration of an upstream lake by reducing the 

supply to the downstream lake from the subwatershed containing 
the upstream lake by multiplying by (1-R!) where R! is the 
phosphorus retention coefficient of the upstream lake. R! 
is calculated as in Step 8, using 9. for the upstream lake. 
The phosphorus load from precipitation, Lop, is 75 mg m=¢ yr! 


Calculate the supply from precipitation as: 
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STEP 14. 


STEP Wd. 


See Ga 


ans ae (kg yr 


If Jy > Joerm’ ie. if the natural supply is greater than 
the permissible supply, allow no (further) development. 
Determine the present number of cottages (N.) and permanent 
dwellings (Ny) within 300 m (1000 ft.) of the lake or any 

of the inflowing streams or rivers from recent aerial 
photographs or field surveys. For cottages, assume 253 
Capita-days per year (0.69 capita-years yr!) in the Kawarthas, 
230 capita-days per year (0.63 capita-years yr!) elsewhere. 


Assume 4.3 people per dwelling, and calculate Ney the number 


of capita-years yr! spent at the lake: 


Ney =e ON63 -} Ne ame ae ee | Ny 


(0.69 replaces 0.63 in the Kawarthas) 

ie. one permanent unit = 6.8 cottage units (6.2 in Kawarthas ) 
Calculate the phosphorus supplied to the lake from the cottage 
units (artificial supply) as: 


] 


Jn = 0.8X Noy (1-R.) (kg yr’) 


A 


where R. = 0 for conventional septic tank - tile field systems 
on the Precambrian Shield. If there is firm evidence that 
holding tanks are used for all household wastes and the systems 
are pumped and removed to a treatment plant outside of the 


watershed, neglect such cottage(s) in the calculations. If 
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there is firm evidence of properly functioning iron-clay 
tile beds, use a figure of 0.16 kg/capita-year instead of 
0.8 kg/capita-year (R = 0.8) for the appropriate locations. 
The artificial load (Ly) is 


3, x 10° 


A 
-2.-1 
A (mg m~ yr’) 
0 


STEP 16. Calculate the present total supply of phosphorus to the lake: 


-] 
Jr = Jy + Jy CKGYE?) 


If J J » allow no further development. 


te perm 


Seale sane k is Jy < Joerm? calculate the total permissable number of 


cottage units: 


ha - J ko/yr 
Noerm = “perm N 
0.63 X 0.8 COUT a ie cane 
OF : cottage yr capita yr 


STEP 18. The additional number of cottage units permitted is: 


Napp ~ Noerm 7 Ney 


“STEP 19. Check to be sure that Noerm does not exceed that BOAT/LIMIT 
figure given by the Lake Alert method. 

STEP 20. With the aid of personnel of the Ministry of Natural Resources, 
locate the cottage units to avoid damage to spawning grounds, 
wildlife habitat, etc. With the assistance of personnel of the 
Ministry of the Environment, locate the septic tile beds in the 
most suitable locations. (the points mentioned above may, at 


the planners discretion, serve as constraints reducing the 


permissible number of cottage units below Noy) 
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EXAMPLES 


EXAMPLE A. To determine the development capacity of SKELETON LAKE, MUSKOKA. 


SEP 
1. Skeleton Lake is ideal for a cold water fishery and body contact 


recreation and should be kept at Level 1 if possible. Therefore: 


[ch] a] = 2 mg n> 


2. Permissible [P] = 9.9 mg fi”? 


3. A = 20°x210°en 


7 2820.0 if 
= 00 ton 
4. Ay = 58) 10° mé 


] 


6. yy © 1245%cRsymi-t="0.500 m yr” 


O2 - Since Ag <10 Ay , determine 


Pr = 38" = 0.965 m 
Ev = 24" = 0.610 m 


"<Q = 58x 10° x 0.500 + (20 X 10°) (0.355) = 36.1 X 10° myr 
5 £290 110° + (2ONx, 10°) (Oneas yh = Ge ye 
600 X 10 
Mag = 2m yo) cs Beale) = ese 
2 20 X 10° 


8. R = 0.426 exp (-0.271 X 1.81) + 0.574 exp (-0.00949 X 1.81) 


0.83 


9. Response time 


= (3+5) (0.17)(0.69 
0.060 


5.9 - 9.8 yr 


eles 


959. % 30 X 0,06 
perm. == 


0.17 


105 mg myn” | 
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4 6 
‘ipempealinie: 20 X10 
10 
= 2100 kg yr”! 


All of Skeleton Lake's watershed is of one geology and land use - 


forested igneous. 


Zw =| 


". £ =4o7 mom yr 


“ de 4.7 X 58 X 10° mg yr"! 
= 273 kg yr! 


= 75 X 20 x 10° 


-] 
kg yr 
10° 


PR 


] 


1500 kg yr_ 


Jy = 273 + 1500 kg yr’ 


1773 kg P yr7! 


N * Joerm 


N. = 370 cottages 
N, = 0 


370 cottage units 


J, = 0.8 X 370 X 0.63 kg P yr! 


A 
= 186.5kg P yr’! 
J, = 1773 + 186.5 kg p yr"! 
= 1960 kg P yr” 
Jr + Joerm 
N = 2100 - 1773 cottage units 


Pern EO Goan 


= 649 cottage units 
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The additional number of cottage units allowed is: 


Napp = 649 - 370 


279 cottage units 


ie. 279 cottages or 41 permanent dwellings 


ALTERNATE CALCULATION 
If it were decided to allow Skeleton Lake to deteriorate to Level 2 


then the following calculations apply: 
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TG. 
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[P] = 18.5 mg m7? 
as before 
Loerm = 18.5 X 30 X 0.060 
Ol? 
5 -] 
Joerm = 3918 kg yr 
-] 
Jp = 273 kg yr 
Jnc = 1500 kg yr! 
PR 
b -| 
Jy = 1/73 kg yr 
Nou = 370 cottage units 
-] 
4 = 186.5 kg yr 
A 
. -] 
Jr = 1960 kg yr 
N = 3918 - 1773 
perm eM a Ee 
OsGs 80 o 


4260 cottage units 


Napp = 3890 cottage units 


3890 cottages or 570 permanent dwellings 
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je. a small town would be permissible 


EXAMPLE B. To determine the effect of a development of 250 cottages on the 
water quality of BOB LAKE, HALIBURTON. 


STEP 
2 20” ie 
me 18.0 m 
Vi. =e 40, FIC Gat 
4. A, = 29.0% 10° wé 


Bie= lees cfs/mi = -43lamyr 
G1 OP = -2OTONe 10> KO 


= 12.5 X 10° myr (A, >10A,) 
o = 0.31 yr’! 

' pam = 

7 Game 25 00K 100 = eS o ye 
2.27 X 10° 
Bek) = 064 
Oe On 26s XA ae ory 
% a0 ona 


Response time = 2.4 - 4.0 yr 


Tig: Coe 7K 2 10° = kg P yr"! 
10 
= 136.3 kg P yr. 
6 
12). Ge 2 ew 2eezeaalO 
PR aa 
10 
= 170.3 kg P yr” 
16. Jy oF Migepse ions 
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307 kg P yr. 
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Ly = 135 mg m™e yr”! 


14. At present 


Ne = 25 cottages 


No = 3 permanent dwellings 
aN Su ee Ome”, ae Es 
bY “Ota. 


45.5 cottage units 


15. J, = 0.8 X 45.5 X 0.63 
= 22.9 Kg P yr" 
16. Jp = 330 kg P yr! 


Oa ie 
L, = 145 mg Pm~ yr 


The predicted spring phosphorus concentration at present is: 


eau (| ah 
Zo 


145 (0.36) 
13. 0¢h0.31) 


= 9.4 mg m7? 
The measured spring total phosphorus concentration (1972) was 


8.5 mg ‘oe very close. 


Predicted summer average chlorophyll] a concentration is: 1.87 mg n> 


(Level 1) cf. measured 1.23 mg nm”? 


Effects of a development of 250 cottages: 


14. Noy = 45.5 + 250 
= 295.5 cottage units 

15. J, = 0.8 X 285.5 X 0.63 
149 kg P yr 
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+ SG - 
L; = 201 mg m™¢ yr! 
The predicted spring phosphorus concentration at the new 


steady-state level is: 


[P] 201 (0.36) 
18.0 (0.31) 


13.0 mg n> 


The predicted summer average chlorophyll] a concentration TS 
3.0 mg ie (Level 2). The Secchi disc readings would average 


about 3 m compared to the present 5 -7 m (see Figure 11). 


EXAMPLE C. To determine the development capacity of MISKWABI LAKE, HALIBURTON. 
STEP 
1. Miskwabi Lake is a deep lake of high water quality and should not 
be developed beyond Level 1. 


[chl a] = 21mg m> 


2. Permissible [P] = 9.9 mg n> 


3c AL = 5 2968 | 10° mé 
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WW = 51.5 X 10° of 
4. Ay = 10.5 X 10° mé 

Uae ay 
5. n =. le22 cfs/mi- = 0.421 myr 
6 Ag < 10 Ao 
- gq = 10.5 xX 10° x 0.421 + 2.68 (0.838 - 0.610) x 10° 
=. 5.03 X10" m 
5 = 5.03% 10° = 0.098 yr 
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82 
(0.18) In 2 
0.098 
1.27 
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= 9.9 X 19.1 X 0.098 
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49.4 kg P yr 
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49 + 201 
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Therefore, allow no further development to maintain the lake at 


Level 1. 


Calculate how many cottage units would result in deterioration 


to Level 2. 
a -3 
[ch] a] = 5 mgm 
id -3 
ml Pitee, 18.57 imgsm 
LoeRM = oc Ona Hon ke OL UOS 
Ost8 
= 192 mg m yr”! 
x -| 
Joepm = Sie Kg Plye 
a -| 
Jy = 250 kg P yr 
; Jn = 266 kg P yr”! is allowed 
ND ERM eee ee 
0.63 X 0.8 
= 528 cottage units (total) 
or an additional 465 cottage units (or 68 permanent dwellings) 


What is the expected steady-state concentration of phosphorus in 


Miskwabi Lake based on the present supply of phosphorus? 
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105 mg m7? yr! 
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10.1 mg m_ 


This compares very favourably to the measured concentration of 


11 mg m73, 
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